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We introduce a novel concept of generalized synchronization, able to encompass the setting of 
collective synchronized behavior for mutually coupled systems and networking systems featuring 
complex topologies in their connections. The onset of the synchronous regime is confirmed by the 
dependence of the system's Lyapunov exponents on the coupling parameter. The presence of a 
generalized synchronization regime is verified by means of the nearest neighbor method. 
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INTRODUCTION 

Synchronization of chaotic systems is a subject widely 
studied in recent years, having both theoretical and ap- 
plied significance [l^. Of a particular interest is the in- 
tricate phenomenon of generalized synchronization (GS), 
originally introduced as an emerging collective motion of 
unidirectionally coupled chaotic oscillators [H, GS has 
been, indeed, observed in numerous systems, both nu- 
merically and experimentally [7H9|, and its main 
features [iQj have suggested many possible applica- 
tions |lll4l3|. 

The very same concept of GS has been introduced ini- 
tially in the case of two unidirectionally coupled oscilla- 
tors, the drive (or master) x(t) system and the response 
(or slave) \i(t) one. GS means the presence of a time 
independent functional relation F[-] between the master 
and slave system states, after a suitable transient time 
interval is elapsed i-e. 



u(t)=F[x(t)]. 



(1) 



In this Paper, we report on the extension of the concept 
of GS for both oscillators with a bidirectional coupling, 
and for complex networks' architectures. It is, indeed, 
evident that Eq. ([1]) needs to be extended for reflect- 
ing the fact that a mutual interaction exists between the 
systems under study. The traditional definition of GS in 
the form of ([T]) introduced for the unidirectionally cou- 
pled chaotic oscillators is based on the fact that the drive 
system state x(t) does not depend on the state of the re- 
sponse system \i(t). Hence, for the GS regime Eq. ([1]) 
may be applied at every moment of time t. In the case of 
the mutual type of coupling between systems the situa- 
tion is radically different, since the interacting oscillators 
are equivalent from the point of view of the coupling, and 
the evolution of the first system is determined not only by 
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the vector x(t), but also by the the state of the second 
system u(t), and vice versa. Therefore, in the general 
case the functional relation between states of the mutu- 
ally coupled systems can not be written in the explicit 
form ([1]), and, as a consequence, the implicit functional 
relation must be used. To take into account the common 
action of the systems on each other, we here propose to 
modify Eq. ([T]) for two chaotic oscillators as 



F[x(t),u(i)] =0. 



(2) 



Note, Eq. (|2]) may be used to inspect both the unidirec- 
tional and the bidirectional types of the coupling (as well 
as any coupling with arbitrary asymmetry between the 
systems). Furthermore, Eq. ([1]) can be considered as a 
special case of (|2]). 

For a generic ensemble of N elements Eq. (j3]) can be 
rewritten as following 



F[xi(t),X2(t),...,x,(t),...,x^(t)] =0, 



(3) 
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where Xi(t) is the vector state of the i-ih element of the 
network. In other words, we assume that the generalized 
synchronization regime implies the presence of a func- 
tional relation between the system states, as before, but 
this relation takes the implicit form (|3]) instead of ([T]), 
i.e. it implies the arousal of a generic manifold in phase- 
space wherein the overall trajectory is lying during its 
time evolution. 



I. GENERALIZED SYNCHRONIZATION, 
LYAPUNOV EXPONENTS AND THE NEAREST 
NEIGHBOR METHOD 

For the unidirectionally coupled oscillators the Lya- 
punov exponents calculation is known to allow to de- 
tect the GS boundary more precisely in comparison with 
the nearest neighbor method, whereas the last one gives 
mostly only the qualitative confirmation of the GS regime 
presence. Therefore, we propose to use the Lyapunov 
exponents to detect the GS regime in the mutually cou- 
pled systems. As far as the nearest neighbor method is 
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concerned, we use it for the verification of the obtained 
results. 

If the dimensions of the drive and response systems are 
Nd and Nr respectively, the behavior of the unidirection- 
ally coupled oscillators is characterized by the Lyapunov 
exponent spectrum Ai > A2 > • • • > Aat^+at^. Due to the 
independence of the drive system dynamics on the be- 
havior of the response one, the Lyapunov exponent spec- 
trum may be divided into two parts [HI, HH : LEs of the 
drive system Xf > ■ ■ ■ > A^^ and LEs of the response one 
Ai > • • • > A^ . For the GS detection in the unidirection- 
ally coupled chaotic oscillators the Lyapunov exponents 
calculated for the response system play the key role (since 
the behavior of the response depends on the drive, these 
Lyapunov exponents are called conditional). With the 
increase of the coupling strength the largest conditional 
Lyapunov exponent A^ becomes negative at the onset of 
the GS regime [l7|. Thus, the negativity of the largest 
conditional Lyapunov exponent 

X[ < (4) 

is considered as a criterion of the GS presence in the 
unidirectionally coupled dynamical systems [sl, [l6| . 

Since in the case of the mutual coupling the spectrum 
of Lyapunov exponents cannot be divided into two parts 
corresponding to the drive and response systems, condi- 
tion dU must also be modified. Notice that, here, the 
dimension of each element is assumed to be = 3, but 
this analytical study may be extended easily to the other 
systems with arbitrary dimensions N^. 

As it has already been mentioned, in the neighborhood 
of any moment of time t [35| Eq. ([2]) may be consid- 
ered as the definition of the implicit functional relation 
between system states, and, therefore, according to the 
implicit-function theorem [isj, locally, the implicit defi- 
nition of the functional relation between system states 
may be used, i.e., x(t) = F[u(t)] or u{t) = F[x(t)]. 
Let us assume without the lack of the generality that 
for — S < t < -\- S (where S is infinitely small) the 
implicit functional relation 

x(r) = F[u(r)] (5) 

is defined with the help of Eq. (j2j). In this case, locally, 
in the given range — S < t < + J, we deal with 
the already studied case of Eq. ([1]). Eq. ([5j) means that 
(under assumption = 3 made above) the following 
local Lyapunov exponents characterize the dynamics of 
the systems: A^ > 0, A^ = 0, A^ < 0, Af 2,3 < 0. In 
other words, in the selected area of the 6D phase-space 
the manifold corresponding to the GS regime is charac- 
terized by one unstable direction and one direction 
with the neutral stability lying inside this manifold, 
whereas all other directions are stable. These directions 
correspond to one positive, one zero and four negative 
Lyapunov exponents. The same statement is also correct 
for the other moments of time t*, although the implicit 
functional relation may take form u(t*) = F[x(t*)] in- 
stead of dSJ. So, one can come to conclusion that in the 



case of the mutual type of coupling the manifold cor- 
responding to the GS regime at every moment of time 
is characterized by one unstable direction, one direction 
with the neutral stability and four stable directions. So, 
having calculated the spectrum of Lyapunov exponents 
for two bidirectionally coupled chaotic oscillators one ob- 
tain that the GS regime in this case produces one pos- 
itive, one zero and four negative Lyapunov exponents. 
Since in the case of the mutual coupling one cannot pick 
out the conditional Lyapunov exponents, the criterion of 
the GS regime should be written in the form 

A3 < 0, (6) 

whereas Ai > 0, A2 = 0. 

Having extended this theoretical consideration to the 
complex networks (where the implicit form of the func- 
tional relation is given by (|3])) one can obtain that condi- 
tion (j6j) remains also correct for the GS regime existence 
in these networks. 

To validate the presence of the generalized synchro- 
nization regime, in parallel with the calculation of the 
spectrum of Lyapunov exponents, one can also make use 
of the nearest neighbor method [ij^lill. The main idea of 
this technique consists in the fact that the presence of the 
functional relation between the interacting system states 
means that all close states ("origins") in the phase sub- 
space of a given system x(t) should correspond to close 
states ( "images" ) in the phase sub-space of the other one 
u{t) (see [14] for details). For mutually coupled oscilla- 
tors the inverse statement must be also correct, i.e. all 
close states in the phase sub-space of the second system 
u{t) must correspond to close states of the first one x(t). 

As a numerical indicator of the existence of a functional 
relationship between the interacting systems, the mean 
distance between images u^'^^ of nearest neighbors x^'^^ 
normalized by the average distance S of randomly chosen 
states of the first system, i.e. 

M-l 

^=M^ 

k=0 

can be calculated (here M is the number of the 
points chosen randomly) [l9|. This characteristic al- 
lows to reveal the qualitative changes in the syn- 
chonous/asynchonous behavior of the coupled systems. 
When the coupling between systems is very small and 
oscillators show the asynchronous dynamics the value of 
this measure is d ~ 1. The phase synchronization regime 
causes the sharp decrease of the value d but it differs from 
zero sufficiently, as before, whereas deep inside the gen- 
eralized synchronization region d tends to be zero due to 
the presence of the functional relation between states of 
the interacting systems. In the vicinity of the GS bound- 
ary (both inside and outside the GS area) the value of 
d-quantity decreases slowly to zero value, with the onset 
of the generalized synchronization regime corresponding 
approximately to the middle of this transition interval. 
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Unfortunately, the nearest neighbor method does not al- 
low to detect precisely the boundary points of the GS 
regime, but it allows to confirm the presence of GS and 
to estimate the location of the GS boundary. 

To illustrate the proposed approach, we use several 
cases, and start our considerations from two mutually 
coupled Rossler oscillators 



^1,2 = -^1,2^1,2 - ^1,2 + £{X2,1 

yi,2 = ^1,2^1,2 + ayi,2, 

Zl,2 =P + ^1,2(^1,2 -C), 



^1,2), 



(8) 



where xi^2(^) = (^1,2,^1,2,^1,2)^ are the vector-states of 
the interacting systems, e is a coupling parameter. The 
control parameter values have been selected by analogy 
with our previous works (2Q|-[22| as a = 0.15, p = 0.2, 
c = 10. The parameter uji^2 defines the natural frequency 
of oscillations, with uji = 0.99 and UJ2 = 0.95 being fixed. 

The behavior of the four largest Lyapunov exponents 
for the considered case is shown in Fig. (TJ a and the de- 
pendence of the quantitative measure ([7]) on the cou- 
pling parameter strength e is given in Fig. [TJ^. One can 
see that at sle ^ 0.106 the second Lyapunov exponent 
A2 passes through zero and becomes negative. There- 
fore, the generalized synchronization regime is expected 
to be observed above the critical value €le- As far as the 
quantitative measure ([7]) is concerned, the curve d{e) de- 
creases monotonically when the coupling parameter value 
increases. The (e; (i)-plane can be divided into 4 parts: I 

- £ G [0; 0.04), the ^-characteristic decreases very sharply 
indicating the transition from the asynchronous motion 
to the phase synchronization regime at eps = 0.04; II 

- £ G [0.04; 0.09), d is not practically changed being 
the evidence of the phase synchronized motion; III - 
e G [0.09; 0.12), the ^-characteristic decreases slowly in- 
dicating the occurrence of the GS regime; IV - £ > 0.12, 
d ^ 0. Note, the bifurcation point €le ^ 0.106 (where 
the second Lyapunov exponent A2 crosses the zero value 
and becomes negative) corresponds exactly to the mid- 
dle of the transition interval III where ^-characteristic 
decreases slowly to zero value. Later, dit els ^ 0.169 the 
lag synchronization (LS) regime takes place. 

In Fig. [Hc-j the phase portraits of the interacting 
Rossler systems ([8]) are shown for the different values 
of the coupling strength s. In the phase portraits of the 
first system x(t) (Fig. [TJ c,e,g,i) three randomly chosen 
points with its nearest neighbors x^"^ are shown. ¥\g- 
xn:Q^^d,f,h,j illustrate the corresponding states u^'^^ in 
the phase sub-space of the second system \x{t). 

One can see easily that for the small values of the cou- 
pling parameter {e = 0.01) all points of the second system 
are distributed randomly throughout the whole attrac- 
tor (Fig. [U d). When the coupling parameter value in- 
creases the points become to be concentrated in a limited 
portion of attractor, with the radius of the distribution 
area decreasing (compare Fig. ^f,h). For e > Ele all 
states of the second system corresponding to the near- 
est states of the first oscillator are also nearest, and vice 



X 


-0.1 
-0.2 
-0.3 

-0.4 
-0.5 

d 
0.8 
0.6 
0.4 
0.2 

























0.05 0.1 ^ 0.15 0.2 s 



I 


























III 

^LE 

-4^ 






















Els 




II 







5/ 

-5' 

-15 



0.05 0.1 ^ 0.15 0.2 s 




FIG. 1: (Color online) (a) The dependencies of the four 
largest Lyapunov exponents on the coupling strength e for 
dH). (b) The quantitative measure d Q versus the coupling 
parameter strength e. The critical values of the coupling pa- 
rameter Sle = 0.106 and sls = 0.169 are marked by arrows. 
(c-j) The phase portraits of Rossler oscillators for the dif- 
ferent values of the coupling parameter: (c-d) e = 0.01 (the 
asynchronous state); (e-/) s = 0.05 (the PS regime); (g-h) 
£ = 0.12 (the GS regime); (i-j) £ = 0.18 (the LS regime). 
Figures {c,e,g,i) show the chaotic attractor of the first sys- 
tem x(t) with three randomly chosen points and its nearest 
neighbors x'^^. Figures {d,f,h,j) illustrate the corresponding 
states u'^''^^ in the phase sub-space of the second system u(t). 
(k) The dependence of the mean distance A between shifted 
states of the interacted systems jS]) 



versa (Fig. ^ g,h and Fig. [TJ ij), that proves the occur- 
rence of GS. Note also the difference between GS and 
LS which consists in the fact that in the LS regime the 
representation points corresponding to the nearest neigh- 
bors are practically in the same part of chaotic attractor 
(Fig.[Tl ij) whereas in the GS regime they can be located 
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in the slightly different regions (Fig. ^ g,h). The addi- 
tional evidence of the fact that the GS and LS regimes 
differ from each other is the behavior of the mean dis- 
tance A between the interacted system states shifted in 
time on the coupling parameter value. Such dependence 
is shown in Fig. [TJA:. It is clearly seen that in the LS 
regime A ^ whereas in the GS regime it is a positive 
one. 



II. GS IN MUTUALLY COUPLED SPATIALLY 
EXTENDED SYSTEMS 

The very same results have also been obtained for 
two mutually coupled Pierce diodes being the classical 
models of beam-plasma systems, demonstrating the com- 
plex spatio-temporal oscillations including the chaotic 
ones [23I, [IJ. The dynamics of Pierce diodes (in the 
fluid electronic approximation) is described by the self- 
consistent system of dimensionless Poisson, continuity 
and motion equations (25| : 



dvio 



(0^1,2) (Pl,2-1), 
S(Pl,2^1,2) 



-^1,2- 



dx 



(9) 



dip 



1,2 



dx 



dt ' dx 

with the boundary conditions 

^1,2(0, t) = 1, pi,2(0,t) = 1, (^1,2(0, t) = 0, (10) 

where (pi^2{x^ t) is the dimensionless potential of the elec- 
tric fleld, pi,2(^,^) and Vi^2{x^t) are the dimensionless 
density and velocity of the electron beam (0 < x < 1), the 
indexes "1" and "2" correspond to the flrst and second 
coupled beam-plasma systems, respectively, ai = 2.8587r, 
a2 = 2.8607r are the control parameters. The bidirec- 
tional coupling between considered Pierce diodes is real- 
ized by the modiflcation of the boundary conditions on 
the right boundar y o f the systems, in the same way as it 
has been done in 



T Q i tne 
|2i,[23 



^1,2(1,^) = e{pi^2{x = l,t) - p2^i{x = l,t)), (11) 

where £ is a dimensionless coupling parameter. 

Continuity and motion equations of (|9|) are integrated 
numerically with the help of the one-step explicit two- 
level scheme with upstream differences and the Poisson 
equation is solved by the method of the error vector prop- 
agation The time and space integration steps have 
been taken as At = 0.003 and Ax = 0.005, respectively. 

For the GS regime detection the nearest neighbor 
method and calculation of the spectrum of Lyapunov ex- 
ponents have also been used. In Fig. [21 a the dependencies 
of four largest Lyapunov exponents on the coupling pa- 
rameter e are shown. For computation of the spectrum of 
spatial Lyapunov exponents the method proposed in [29| 
has been used. It is clearly seen that as in the case of 




mutually coupled Rossler systems two Lyapunov expo- 
nents do not practically depend on the coupling param- 
eter £, i.e. one Lyapunov exponent Ai is always positive 
(except the windows of periodicity) whereas the second 
one A3 remains zero. At the same time, two Lyapunov 
exponents (initially positive A2 and initially zero A4) de- 
pend on the coupling parameter and pass in the fleld of 
the negative values with the coupling parameter value 
increasing. One can assume, that as in the case of mutu- 
ally coupled Rossler systems the transition of the positive 
Lyapunov exponent A2 in the fleld of the negative values 
(for e = Ec = 0.078) is connected with the generalized 
synchronization regime onset in mutually coupled beam- 
plasma systems. 

To conflrm the assumption made above the near- 
est neighbor method has been used. To characterize 
the degree of closeness of the interacted system states 
quantitatively the measure d deflned by ([7j) has also 
been computed. As the interacted system states vec- 
tors ui^2{x^t) = (pi,2, '^1,2, ^^1,2)^ have been used, whose 
norm 1 1 • 1 1 has been calculated as 



(12) 



In Fig. [21 6 the dependence of the quantitative measure d 
on the coupling parameter e is shown. It is clearly seen 
that (^-characteristics decreases monotonically from one 
to zero with e value increasing. At that, Sc is approxi- 
mately in the middle of the falling fleld e G [0.04; 0.12], 
that indicates the generalized synchronization regime 
presence. It should be noted that generalized synchro- 
nization does not coincide with the complete synchro- 
nization regime in this case. For the control parameter 
values mentioned above it is realized for e ~ 0.17. 

Additional evidence of the presence of generalized 
synchronization regime in two mutually coupled beam 
plasma systems is the behavior of the nearest neighbors 
in the phase space of interacted systems. In Fig. ^c-h 
the reconstructed attractors of interacted Pierce diodes 
on {pi^2{x = 0.2,t)^ pi^2{x = 0.6, t))-plane for different 
values of the coupling parameter e are shown. On at- 
tractors of the flrst system {Fig.^c,e,g) three randomly 
chosen points and its nearest neighbors are also indi- 
cated. Fig. ^dj,h illustrates the corresponding states 
in the phase space of the second system. 

It is clearly seen that as in the case of Rossler systems 
for small values of the coupling parameter {e = 0.002) 
all points in the phase space of the second system are 
distributed randomly over all attractor (see Fig. [21 d). 
With the coupling parameter value increasing the points 
begin grouping in the limited range of attractor that 
corresponds to the phase synchronization regime onset, 
with the radius of such fleld being decreased (compare 
Fig. ^f,h). For e > Sle all states of the second beam- 
plasma system corresponding to the nearest neighbors 
of the flrst Pierce diode are also nearest and vise versa 
(Fig. ^g,h)^ that is the evidence of the generalized syn- 
chronization regime presence. 
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FIG. 2: (Color online) (a) The dependencies of the four 
largest Lyapunov exponents on the coupling strength £ for (|9|) . 
(b) The quantitative measure d ([7|) versus the coupling param- 
eter strength s. The critical value of the coupling parameter 
£le = 0.078 corresponding to the zero-cross of the positive 
Lyapunov exponent is marked by arrow, (c-h) The recon- 
structed attractors of two mutually coupled Pierce diodes on 
(pi,2{x = 0.2, t), pi,2{x — 0.6, t))-plane for the different values 
of the coupling parameter: {c,d) e = 0.002 (the asynchronous 
state), (ej) £ = 0.05 (the PS regime), {g,h) £ = 0.10 (the 
GS regime). Figures (c,e,g) show the reconstructed attractor 
of the first Pierce diode with three randomly chosen points 
and its nearest neighbors. Figures (d,f,h) illustrate the corre- 
sponding states in the phase sub-space of the second Pierce 
diode 
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FIG. 3: (Color online) The dependencies of the four largest 
Lyapunov exponents on the coupling strength £ 



So, on the basis of consideration carried out one can 
conclude that in both cases considered above, the onset 
of the GS regime is connected with the sign change of the 
initially positive Lyapunov exponent A2 in the same way 
as for two unidirectional oscillators \r?\. As well as in the 
case of the chaotic oscillators coupled unidirectionally, 
for the systems with the bidirectional coupling the onset 
of the GS regime €le precedes the boundary of lag or 
complete synchronization sls- 



III. GS IN MUTUALLY COUPLED LORENZ 
OSCILLATORS 

To prove the made decisions conclusively, in this Sec- 
tion of the paper we consider another example of the os- 
cillators coupled mutually, namely, two Lorenz systems 



^1,2 = cr(^l,2 - 

yi,2 = 22:1, 2 

^1,2 = -bzio -{ 



^1,2) + ^(3:2,1 - ^1,2), 

- yi,2 - ^1,2^1,2, 
^1,2^1,2. 



(13) 



where xi^2(^) = (^1,2,^/1,2,^1,2)"^ are the vector-states 
of the interacting systems, £ is a coupling parameter, 
a = 10.0, b = 8/3, ri = 40.0 and r2 = 35.0. Due to the 
bistable type of the chaotic attractor of the Lorenz os- 
cillator there are certain particularities of the GS regime 
which do not take place in the systems with the Rossler- 
like chaotic attractor. 

According to the decision made in Sec. HI HH the onset 
of the GS regime is connected with the sign change of 
the initially positive Lyapunov exponent A2 at e^^; ^ 6.0 
(Fig. [3|). To prove this statement for two mutually cou- 
pled Lorenz systems ([13]) we have used the nearest neigh- 
bor method again both below (Fig. [H a,b) and above 
(Fig. m c,d) the critical point €le- For this purpose the 
reference point and its nearest neighbors x^"^ have 
been selected in the phase space of the first Lorenz system 
(Fig. m a,c) and the corresponding to them points u^'^^ 
have been found in the the phase space of the second 
Lorenz oscillator (Fig. [3J b,d). For e > Sle the existence 
of the GS regime is evidenced by the fact that all states 
of the second Lorenz system corresponding to the near- 
est states of the first oscillator are also nearest (Fig.[4l d). 
Alternatively, below the threshold e the points of the sec- 
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FIG. 4: (Color online) The phase portraits of two mutually 
coupled Lorenz oscillators fT3)) for e — 5.7 (a^^) and e = 
6.1 (c,d). Figures {a,c) show the chaotic attractor of the 
first system x(t) with the reference point and its nearest 
neighbors x^^. Figures {b,d) illustrate the corresponding to 
them states the phase space of the second system 

u{t) 



ond system are located on the both sheets of the chaotic 
attractor (Fig. IH b) that indicates the break of GS. 

Again, as well as for two mutually coupled Rossler sys- 
tems (jSj) and Pierce diodes (j9j), the onset of the GS regime 
is shown to be connected with the sign change of the ini- 
tially positive Lyapunov exponent A2 taking place when 
the coupHng strength e grows. 

At the same time, in the dynamics of mutually cou- 
pled Lorenz oscillators (p!3|) there are certain particular- 
ities connected with the occurrence of GS caused by the 
bistable type of the chaotic attractor of the system under 
study [36| . One can see easily that in the case under con- 
sideration the coupling strength value corresponding to 
the onset of GS is larger sufficiently than the analogous 
value for two mutually coupled Rossler systems ([8]) and 
Pierce diodes (j9j). Owing to the great coupling strength 
value, in the vicinity of the GS onset (below the bifur- 
cation point Ele) two interacting Lorenz systems are 
greatly synchronized with each other almost all time ex- 
cept for the short time intervals when the representa- 
tion point of one of the coupled oscillators remains in 
the one sheet of the chaotic attractor whereas the repre- 
sentation point of the second oscillator jumps to another 
sheet (see Fig. [5]). After such a short jump both phase 
trajectories approach each other, and the oscillators start 
showing synchronous dynamics again. It is the short-time 
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FIG. 5: (Color onhne) (a) The fragment of the time series of 
two mutually coupled Lorenz oscillators (p^ corresponding 
to the short-term period of the phase trajectory divergence. 
(6) The phase portraits of Lorenz systems and the phase tra- 
jectories corresponding to the time interval shown in Fig.[5l a. 
The coupling strength s = 5.7, the GS regime is not observed 
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FIG. 6: (Color online) (a) The fragment of the time series of 
two mutually coupled Lorenz oscillators ()13p corresponding 
to the short-term period of the phase trajectory divergence. 
(6) The phase portraits of Lorenz systems and the phase tra- 
jectories corresponding to the time interval shown in Fig. [6l a 
The coupling strength e — 6.1, the GS regime is detected 



phase trajectory divergence that is responsible for the GS 
regime does not come into being. Above the critical point 
£le there are also the phase trajectory divergencies, but 
they do not envelop both sheets of attractors and the 
representation points remain in the limits of one and the 
same sheet during this perturbations (Fig. [6|). 
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The difference between these types of dynamics ob- 
served below and above the critical point Sle allows to 
explain the occurrence of the GS regime in two mutu- 
ally coupled Lorenz systems (fT3]) . When the representa- 
tion points are located at one and the same sheet of the 
chaotic attr actor, the functional relation between vec- 
tor states ([2]) is likely to exist, since the observed types 
of behavior is very close to the complete synchroniza- 
tion regime due to the rather large value of the coupling 
strength e. Alternatively, the divergence of the phase tra- 
jectories on the different sheets of the chaotic attractor 
terminates this functional relation, and, in turn, the gen- 
eralized synchronization regime is destroyed. So, below 
the critical point sle the GS regime does not take place 
due to the presence of the short-term time intervals with 
the divergence of the phase trajectories on the different 
chaotic attractor sheets. More specifically, the intermit- 
tent behavior near the onset of the GS regime is observed 
(as well as in the vicinity of the other t ype s of chaotic syn- 
chronization like lag synchronization |3Q| , Isif , phase syn- 
chronization [32I, [33I and for generalized s ync hronization 
in the unidirectionally coupled oscillators [3^ ) , that may 
be considered as the additional evidence of the correct- 
ness of the obtained results. Note, that above the onset 
of the GS regime (sle) there are also the phase trajec- 
tory divergencies (which do not destroy the GS regime, 
since the representation points remain in the limits of 
one and the same attractor sheet and, as a consequence, 
the functional relation ([2]) takes place) preventing the oc- 
currence of LS. These perturbations of the synchronous 
dynamics vanish above the critical point sls where the 
lag synchronization regime comes into being. 

So, having considered the behavior of two mutually 
coupled Lorenz system, we have obtained the additional 
evidence of the correctness of the proposed viewpoint on 
the GS regime in the systems with the mutual type of 
the coupling. 

IV. GS IN NETWORKS OF COUPLED 
OSCILLATORS 

Now we move to a more complicated situation, and 
analyze the GS in complex networks. As we have men- 
tioned above, in this case between the interacting system 
states the functional relation in the form (j3j) should be 
established. 

Developing the concept of GS in the mutually cou- 
pled oscillators for the networks, one can say that the 
phenomenon of GS in the complex network can be un- 
derstood as the state of the whole network when the 
co-ordinates of all oscillators consisting this network are 
uniquely determined by the values of co-ordinates of only 
one node x/^ (chosen arbitrary). Following the arguments 
given in Sec. HI one can use for the generalized synchro- 
nization regime the implicit form of the functional rela- 
tion between network's node states 

x,(r) = F[xfc(r)], yij^k, (14) 



where t*— 5<t<t*-\-5 (where 5 is infinitely small). 
Again, locally, in the given range — (5 < t < + (5, 
we deal with the already known case. Under assumption 
Nd = 3 made above the following local Lyapunov expo- 
nents characterize the dynamics of the systems: > 0, 
= 0, A3 < 0, A*i 2,3 < 7^ /c. In other words, in 
the selected area of the 3A^D phase-space the manifold 
corresponding to the GS regime is characterized by one 
unstable direction and one direction with the neutral 
stability lying inside this manifold, whereas all other 
directions are stable. These directions correspond to one 
positive, one zero and (3A^ — 2) negative Lyapunov expo- 
nents. The same statement is also correct for the other 
moments of time t*, therefore, for the complex network 
the manifold corresponding to the GS regime at every 
moment of time is characterized by one unstable direc- 
tion, one direction with the neutral stability and {3N — 2) 
stable directions. So, having calculated the spectrum of 
Lyapunov exponents for the network of coupled chaotic 
oscillators one obtain that the GS regime in this case 
produces one positive, one zero and (37V — 2) negative 
Lyapunov exponents. 

To prove the theoretical assumptions mentioned above 
we consider a network consisting of = 5 Rossler sys- 
tems with slightly mismatched cj-parameter values. The 
evolution of i-th node(z = 1, . . . , A^) is described by the 
following equations 

yi = ujiXi ^ ayi, (15) 
Zi =p^ Zi{xi -c), 

where the values of the control parameters a, c have 
been chosen to be the same of the case of two cou- 
pled oscillators (|8]), uJi = 0.95, UJ2 = 0.9525, UJ3 = 0.955, 
UJ4 = 0.9575, CJ5 = 0.96, Xi(t) = (x^, Zi)^ is the vector- 
state of the i-th node, e is the coupling strength between 
nodes, Gij is the element of the coupling matrix G. G 
is a symmetric zero row sum matrix, with Gij {i 7^ j) 
being equal to 1 whenever node i is connected with node 
j and otherwise, and Gu = — Xlj/i ^ij- The topology 
of the links between nodes in the network under study 
has been selected in such a way that each element of the 
network is connected with each other. 

The dynamics of the considered network is character- 
ized by 3 A" = 15 Lyapunov exponents. If the coupling 
between nodes is equal to zero, there are A" positive, 
A' negative, and A" zero Lyapunov exponents. With 
the increase of the coupling strength e the zero Lya- 
punov exponents as well as the positive ones go grad- 
ually to the region of the negative values. The depen- 
dencies of the seven largest Lyapunov exponents on the 
coupling strength e for the network consisting of five 
Rossler systems are shown in Fig. [71 One can see that at 
Sle ^ 0.0385 the second Lyapunov exponent A2 passes 
through zero and becomes negative. Therefore, the gen- 
eralized synchronization regime is expected to be ob- 
served above the critical value £le- 

To prove the presence of GS we have used the nearest 
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FIG. 7: (Color online) The dependencies of the seven largest 
Lyapunov exponents on the coupling strength s for the net- 
work consisting of five Rossler systems (|15p . The onset of the 
GS regime in the network is shown by the arrow 

neighbor method in the same way as for two mutually 
coupled Rossler systems. In Fig. [8] the phase portraits 
of all Rossler systems of the network are shown for two 
values of the coupling strength, below (Fig. [51 a, £ = 0.03) 
and above (Fig. [HI e = 0.04) the critical point ele- 

In the phase portraits of the three systems Xi(t), 
i = 2 ^ 4 three points (one point for each system) with 
its nearest neighbors have been selected randomly (♦ — 
i = 2,H z = 3, — i = 4) and the points correspond- 
ing to them have been detected in all other coupled sys- 
tems. For 6 = 0.03 (Fig. [3 a) the points are concentrated 
in a limited range of attractor and distributed along the 
radius being the evidence of the presence of PS and the 
absence of GS. For e > Ele (Fig. [HI b) all states of all os- 
cillators are nearest neighbors, thus proving the existence 
of GS. 

Similarly to the case of two mutually coupled Rossler 
oscillators considered above we compare the onset of the 
GS and LS in the network ([T5]) . But due to the small 
values of the control parameter detuning the LS regime 
is very close to the complete synchronization (GS) one. 
In the simulations, the onset of the GS regime can be 
monitored by looking at the vanishing of the time average 
(over a window T) synchronization error 

1 

In the present case, we adopt as vector norm ||x|| = 
^/x^^^^y^^^^z^ . Fig. [9] reports the synchronization er- 
ror versus e for a given topology. One can see that the 
synchronization error becomes close to zero considerably 
later the GS regime arising. When the GS regime takes 
place (E) is still positive that is the evidence of the GS 
regime absence. 

So, the GS regime in networks of coupled nonlinear 
elements can be detected by the moment of transition 
of the second (positive) Lyapunov exponent in the field 
of the negative values. Now we analyze the influence of 
the number of elements and topology of the network on 
the GS regime onset. In Fig. [10] the boundaries of the 



GS regime on the "number of elements — coupling 
parameter -plane for networks of different topologies 
of links between nodes are shown. Gurve 1 corresponds 
to the random network whereas curves 2 and 3 refer 
to the regular and "small- world" networks, respectively. 
For all considered cases the values of the control param- 
eters uji have been selected randomly in such a way that 
the probability distribution density of Co'i-values has been 
obeyed by the Gaussian distribution with the mean value 
ujQ = 0.95 and variance Auj = 0.017 that corresponds to 
the case of the relatively large values of the control pa- 
rameter detuning. It is clearly seen from Fig. [10] that 
the topology of the network influences sufficiently on the 
GS regime onset. In particular, the threshold of the GS 
regime onset decreases for the random network, whereas 
both for the regular and "small world" ones it increases 
monotonically. At that, the boundary of GS for regu- 
lar network grows more rapidly in comparison with the 
"small world" one. 



V. MECHANISMS OF GS OCCURRENCE 

Now, following the approach of some of our previous 
works [3], we move to reveal the mechanisms associated 
with the emergence of the GS regime in the case of a 
generic ensemble of coupled systems (with, obviously, the 
case of two coupled oscillators considered as the simplest 
variant of such a configuration). 

For the purpose of exemplification, and with- 
out lack of generality, we here-below character- 
ize the state of the network by the only vector 
JJ = {ui,U2,...,Ui,...,UN.NdV, where usi-2 = Xi, 
usi-i = i/i, usi = Zi^ instead of the set of vectors 

= {xi^yi^ Zi)^ ^ i = 1,A^. The dimension of each el- 
ement of the network is assumed to be = 3 again, but 
this analytical study may be extended easily to the other 
systems with arbitrary dimensions N^. 

Following the above formalism, the entire network may 
be considered as a high-dimensional autonomous dynam- 
ical system, whose evolution equation is given by 

U = L(U)+£GU. (17) 

Here the vector function L(-) determines the evolution of 
the elements of the network in the absence of the cou- 
pling, whereas the additive term eG describes the influ- 
ence of the topology and the coupling strength of the 
links between oscillators. Matrix G specifies the struc- 
ture of the dissipative couplings between nodes, and it 
is assumed to be a symmetric zero row sum matrix, 
Gii = -^j^iGij, with Gij (i ^ j) being equal to 1 
whenever variable Ui forces the variable Uj and other- 
wise. 

It is easy to see that the term sGU brings the addi- 
tional dissipation into the system ([T7|) . Indeed, the phase 
flow contraction is characterized by means of the vector 
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FIG. 8: (Color online) The phase portraits of five Rossler oscillators for two different values of the coupling parameter: (a) 
£ = 0.03 (the PS regime) and (6) e = 0.04 (the GS regime) 
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FIG. 9: {E) vs e for the network of Rossler oscillators ()15p . 
The onset of the GS regime in the network sle '^^ shown by 
the arrow 




also negative and the dissipation in the considered group 
increases with the growth of the coupHng strength £, re- 
sulting in the simplification of the otherwise chaotic dy- 
namics of the system (pT|) . 

To characterize the complexity of the motion, the spec- 
trum of Lyapunov exponents is frequently used. In 
the case under study, let's suppose that the behav- 
ior of the system (pT|) is initially described by the set 
Ai > A2 > • • • > Aatat^ of Lyapunov exponents, with N 
of them being positive (i.e., while no synchronous motion 
is set up, each one of the elements of the network con- 
tributes independently with one positive Lyapunov expo- 
nent to the overall spectrum). As dissipation increases, 
some of the initially positive Lyapunov exponents be- 
come negative, and each passage of a Lyapunov exponent 
trough zero testifies that one more degree of freedom of 
the chaotic motion corresponds to a contractive direction. 
When A2 becomes negative, only one degree of freedom 
is representative of the evolution of the network, i.e. a 
GS regime is built. Indeed, as soon as A2 is negative, 
all systems have to arrange their evolution into a specific 
collective motion, wherein the functional relation (|3]) is 
taking place. 



FIG. 10: (Color online) Boundaries of the GS regime on the 
"number of elements N — coupling parameter e" -plane for 
random (curve i), regular (curve 2) and "small world" (curve 
3) networks 



field divergence 



NdN 



lim lim — r— = div L + £ ^ ^ Gu , 



At^o Ay^o AV At 



(18) 



where /SV is the elementary volume of the phase space of 
the system (pT|) . Since Gu < 0, the term sY^^J^ Gu is 



CONCLUSION 

In conclusion, we have analyzed the GS regime in sys- 
tems with a mutual type of the coupling. We have ex- 
tended the definition of GS, for being valid also for pairs 
of mutually coupled chaotic oscillators, and for the com- 
plex networks. The GS regime onset in 3D systems is 
shown to be connected with the zero-crossing of the sec- 
ond Lyapunov exponent. GS may be, therefore, consid- 
ered in terms of the transition from the high-dimensional 
hyperchaotic regime to the chaotic oscillations. The ob- 
tained results are proved by means of the nearest neigh- 
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bor method. Since the developed theory is apphcable to 
different systems, we expect that the very same mech- 
anism win be observed in many other relevant circum- 
stances. Particularly, the obtained results could be ex- 



tended to the systems which dimension of the phase space 
Nd > 3 including the spatially extended media and cou- 
pled systems with a different dimension of the phase 
space. 
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